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Abstract 

We classify the Hamiltonians H = p^+Py^ V{x, y) of all classical 
superintegrable systems in two dimensional complex Euclidean space 
with second-order constants of the motion. We similarly classify the 
superintegrable Hamiltonians = + Jf + J3 + V{x, y, z) on the 
complex 2-sphere where + y"^ + = 1. This is achieved in all 
generality using properties of the complex Euclidean group and the 
complex orthogonal group. 
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1 Introduction 



It is known from classical mechanics that a mechanical system with n degrees 
of freedom is completely integrable if there are n functionally independent 
constants of the motion which are mutually in involution The idea of a 
superintegrable system is that there exist more than n functionally indepen- 
dent constants of the motion, but not necessarily in involution. If there are 
2n — 1 such constants the system is said to be maximally superintegrable or 
just superintegrable ^ ^. Here we consider only the case where there 
exist 2n — 1 functionally independent constants of the motion (including the 
Hamiltonian) that are quadratic in the momenta. Ranada investigated 
such systems and noted that many could be found in Drach's list of potentials 
admitting constants cubic in the momenta.^ In the papers 0, ^ we have 
given a complete classification of all non- degenerate potentials on complex 
Euclidean 2-space and on the complex 2-sphere that give rise to superinte- 
grable systems. (For example in we have calculated all the inequivalent 
superintegrable potentials V that are non-degenerate in the sense that they 
depend uniquely on four arbitrary parameters, i.e., one can prescribe the 
values of V,Vx,Vy, Vyy arbitrarily at any regular point (xq, yo) and these val- 
ues determine V{x,y) uniquely.) In this article we relax this requirement 
and ask the same question but without the condition of non-degeneracy: For 
which potentials in two dimensions do there exist at least two constants of 
the motion 

Aj = aj{x,y)pl + bj{x,y)pl + Cj{x,y)p-,Py + dj{x,y) = A'j + dj, j = 1,2 (1) 
in addition to the Euclidean space Hamiltonian 



i.e., {H,Aj} = 0, j = 1,2, and such that the 2n — 1 = 3 constants of the 
motion H = ^0,^1,^2 are functionally independent on phase space? We 
will do the same for the Hamiltonian on the complex 2-sphere 



where x'^ + y'^ + z'^ = 1 and Ji = yp^ - zpy, J2 = zp^ - xp^, J3 = xpy - yp^. 
We give a complete solution. The computations are lengthy, and throughout 
we have made use of a computer algebra package. We give many details in 
the first few examples, to make our method clear. 

^These are [E1,E2,E7,E9,E16,E19,E20] in out notation. 



H = pl+pl + V{x,y) 



(2) 




(3) 
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In references 0, || we worked out the structure of the quadratic algebra 
for each of the non-degenerate potentials. In this paper we supplement those 
results by computing the quadratic algebras for the non-constant degenerate 
potentials. Also we correct a few errors and fill in some gaps in those earlier 
papers. 

2 Superintegrability in E2^c 

In the computations to follow, quite often we will be considering systems 
such that the corresponding Hamilton Jacobi equation H = E, (j^), can be 
solved by the method of separation of variables. Then px = ^,Py = and 
there is a complete integral of the form 

S = U{u,E,X) + V{v,E,X) 

for separable coordinates u = u{x,y),v = v{x,y) and some separation con- 
stant A. This will not always be the case, but when separation is possible 
a knowledge of the separable coordinates u, v will greatly simplify our com- 
putations. In general we will use the structure of the complex Euclidean 
group E{2, C) and its Lie algebra e(2, C) to solve our problem. The ele- 
ments Li = L2 = Py and L3 = M = xpy — yp^ form a basis for e(2, C) 
under the Poisson bracket, and quadratic elements LiLj, i < j form a basis 
for all purely quadratic functions A' such that {A',pl + Py} = 0. Thus the 
quadratic integrals Aj can be written in the form 

Aj = afLkLe + dj{x, y) = A'j + dj (4) 

for suitable constants a^^ = a^^. 

We will not regard Euclidean space Hamiltonians as essentially different 
if they are related by a Euclidean transformation. Because of this we can use 
Euclidean group transformations to simplify the expressions for the Aj and 
classify them into equivalence classes. If we do this then there are equivalence 
classes of constants whose typical representatives are [Q 

pI, {Px + iPyf, M^, M{px + ipy) + {px - ipyf, (5) 

M^ + {Px + iPyf, Mp,, M^ + cY,, M{px + ipy). 

(Note that, up to the addition of an arbitrary multiple of the Casimir el- 
ement pI. + P^, this is simply a choice of a representative on each distinct 
orbit of second order elements in the enveloping algebra of e(2, C) under the 
adjoint action of E{2,C).) Without loss of generality, we can assume that 
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A[ coincides with one of these representatives and use the defining relations 
{H, Aj} = for constants of the motion, to determine the general form of the 
second constant. This is a particularly useful strategy since all but the last of 
the list of representatives has a form that implies separation of variables 
in at least one coordinate system. If we use this fact then we see that the 
corresponding potential V must have the form implied by separation, and 
the requirement of an extra constant of the motion implies strong conditions 
on this functional form. This will greatly simplify our computations. For 
all but one potential, we find that the associated constants determine more 
than one separating coordinate system. These are listed in Appendix 

We now deal with each of these cases individually. Consider the first con- 
stant in our complete family of equivalence classes. We assume our Hamil- 
tonian has a constant of the form 

Ai=pl + d^{x,y). (6) 

The condition {H, Ai} = implies dydi{x, y) = 0, dxiV — di{x, y)) = 0, so 
we can assume that the Hamiltonian has the form 

H = Ao=pl+pl + fix) + h{y) (7) 

where di{x,y) = f{x). For superintegrability we must have one additional 
constant of the motion which can be written 

A2 = 4^LkLe + d2ix,y). (8) 

Since we can always add linear combinations of Aq and Ai to A2 without 
changing the system, we can assume that 

A2 = AM^ + BMp^ + CMpy + Dp^py + d2{x, y) (9) 

where A, B, C, D are constants, not all zero. 

Because of the forms of Aq and Ai we can always apply translations to A2 
in order to simplify its form. Suppose A 7^ in (^. Then, normalizing so that 
A = 1, by appropriate translations in x and y we can pass to a new Cartesian 
coordinate system in which B = C = 0. The condition {^0,^2} = then 
determines the possible forms of A2. Indeed, equating coefficients of px and 
Py, we find 

dyd2{x,y) = {D-2xy)f{x)+2x^h'{y), d^x^y) = 2y'f{x) + {D-2xy)h'{y). 
Equating the cross partial derivatives of d2{x,y) we obtain the condition 

(/" + -/') - {h" + -h') = ^{h" - f"). (10) 
X y 2xy 

If D = 0, the variables separate and we find a well known non-degenerate 
superintegrable potential. 



4 



[El] V = u\.^ + y^) + ^^+l 

The additional constant has the form 

This example together with its Poisson bracket relations is well studied 

If y4 = but 5^ + 7^ in (|^) then we can rotate coordinates so that 
5 = 0, normalize so that C = 1, translate to obtain D = 0, and find the 
non-degenerate potential 

[E2] V = uj^{Ax'^ + y^) + ax+ ^ 

y 

The additional constant has the form 

A2 = Mpy + d2ix,y) . 



There are no further non-degenerate potentials separating in Cartesian co- 
ordinates. We now return to the case A ^ and suppose that D ^ 0. Then 
from (0) we see that h", f" must satisfy a functional equation of the form 
h"{y) — f"{x) = xy{G{y) + F{x)) for some functions G, F. Solving this equa- 
tion and substituting back into (0)) "^^ that the variables separate and 
we obtain the solution 

[E3] V = uj\x^ + y^). 

In this case, di{x,y) = uP'x^, and there are two additional constants, 
one of which is first order. They can be taken in the form 

M = PxPy + uj'^xy , X = M . 

The Poisson bracket relations for these constants are 

{Ai, X} = 2A2 , {^2, X} = Ao- 2Ai , {Au A2} = -2uo^X . 

Since Aq, Ai, A2 are functionally independent, all constants of the mo- 
tion are functions of these. It is easy to verify that X satisfies the 
functional relation 

Al - Ai{Ao - Ai) + cj^X^ = 0. 
This example is the harmonic oscillator in two dimensions. 
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If A = and ^2 + = but S ^ 0, we can take S = 1, C = i (by 
mapping y to —y if necessary), and translate to get D = 0. A straightforward 
computation gives 

[E4] V = a{x + iy). 

Here di{x,y) = ax and the Hamiltonian admits two extra constants, 
one of which is first order, 

A2^M(p^ + ipy) + ^-^(x + iyf, X^p^ + ipy. 
The Poisson bracket relations take the form 

{Ai,X} = -a, {A2, X} = {^1,^2} = -iX^+2iAiX -iAoX, 

with the functional relation 

Al + X^{2Ao - 4^1 + X^) + AiaA2 = . 

\iA = B = C = Q then we normalize D = 2 and find the non-degenerate 
potential V = a;^(a;^ + y"^) + ax + Py. For a fixed choice of the parameters, 
the Hamiltonian admits the first order constant of the motion 

X = 2uj^M + apy - (5p^. 

By an appropriate translation we can obtain X' — M, which is case [E3] . 

There are special cases of potentials [El] and [E2] such that the Hamil- 
tonian admits more than two constants of the motion. The possibilities are 

[E5] V^ax. 

Since in this case, Aq — Ai = Py, we can replace Ai with the first order 
constant Py and take the additional constants as 

A2 = Mpy : ^3 = PxPy + 2 ^ ' ^ 

They satisfy the Poisson bracket relations 

{A2,X}^As, {As,X}^--, {A,,A2}^2X'-AoX 

and the functional relation 

Al + X^- AoX^ + aA2^0. 
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[E6] V=-. 

As for the previous case, we can replace Ai with py. The additional 
constants are 

A2 = MP, A3=M2 + ^, X=Py. 

Their Poisson brackets are 

{A2, X} = Ao-X\ {As, X} = 2A2 , {A3, A2} = -2XAs - 2aX , 
and they satisfy the functional relation 

Al - AsiAo - X^) + aX^ = . 

This concludes the list of possible potentials corresponding to the first equiv- 
alence class of second order elements in the enveloping algebra of e(2, C). 
For orbits of the second type, the constant of the motion Ai has the form 

Ai=p^^ + di{x,y). (11) 

(We adopt the notation p± = px± ipy, z = x + iy,z = x — iy.) It follows from 
the relation {Aq, Ai} = that for a Hamiltonian to admit a constant of the 
motion of the form (|llD the potential V must have the form 

V = f{z)z + h{z), 

for some functions / and h. We can assume 

A2 = AM"^ + BMp+ + CMp_ + Dp\ + d2{z, z). 

There are several possibilities. In the first case we assume A 7^ 0. Then we 
can normalize A = 1, translate to get i? = C = 0, and write D = We 
find the non-degenerate potential ||^ 

r , , , az dz 

[E7] V = + _ — + -fzz. 

\/ z^ — (? \' z'^ — {z + V z'^ — c'^y 

Here the second constant of the motion can be taken in the form 

A2 = + c^pl + d2{x,y). 

The limiting case of this as c — > gives the potential |0 
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Here the second constant of the motion has the form 

A2 = + d2{x,y) . 

If A = but BC 7^ 0, we can normahze and rotate to obtain B = —i/2, C 
i/2, and translate to get Z) = 0. We obtain the non-degenerate potential [] 

[E9] l^=4 + /?x+^^" + '^ 



'Z ^JZ 

The second constant of the motion is 

A2 = Mpy + d2{z, z). 

\{ A = BC = but C 7^ 0, and D ^ we normalize C = 4i and rotate so 
that D = 1 to obtain the non-degenerate potential 

[ElO] V = az + p{z- ^f) + -i{zz - ]^z^). 

Here the second constant of the motion has the form 
A2 = AiMp_ + p+ + d2{x, y). 

If A = 5C = but 5 ^ 0, (or if C 7^ 0, = and D = and we reflect 
y —>■ —y), we can normalize B = 1 and translate so that D = to obtain the 
non-degenerate potential 

[Ell] V = az + ^ + ^. 

\/z 



Here the second constant of the motion has the form 

A2 = Mp^^d2[x,y). 

There are special cases of potentials [E7,E8,E11] that admit two extra con- 
stants of the motion. In each of these cases A\ = jPl-, i.e., di{x,y) = and 
hence X = p_ is a constant of the motion. The possibilities are 

[E12] V 



with the constants of motion given by 



X = p_, A2 = M'--pX- ^ _^ , A3 = Mp_ + 
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The Poisson bracket relations are 

{X, A2} = 2iAs , {X, As} = iX^ , {A2, As} = -2iXA2, 
with the functional relation 

Al-X^A^-'^Al+'^^O. 



[E13] V^-^ 



with the constants of motion given by 

The Poisson bracket relations are 

{X, A2} = iAo , {X, As} = iX'' , {A2, As} = -2iXA2 , 
with the functional relation 



AsAo - X^A2 - 4q;^ = . 
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[E14] V^^ 



with the constants of the motion given by 

X = p_ , ^2 = Mp_ - — , As^M^ + — . 

z z 

The Poisson bracket relations are 

{X, A2} = iX^ , {X, As} = 2^2 , {A2, As} = 2iXAs , 
with the corresponding functional relation 

il - As. 



Al - AsX^ + aAn = . 



[E15] V = h{z) 

where h is any function of z, not necessarily as already given above. A 
constant of the motion always exists of the form 

i f dh 



A I f an . 



in addition to the constant X = p_. Indeed we might take h{z) = az^, 
in which case 

A2 = Mp_ + -az^. 

3 

This is an example of a potential for which separation of variables 
occurs in only one coordinate system, 

For orbits of the third type, the constant of the motion Ai has the form 

A, = M^ + diix,y). (12) 

In this case V = f{r) + where r,9 are polar coordinates (see Appendix 
and di{x,y) = h{0). We can assume that the second constant takes the 
general form 

A2 = AMp+ + BMp. + Cpl + Dp^ + d2(x, y). (13) 

If AB 7^ then we can rotate and normalize to get A = —B = —i/2 and 
translate to achieve C = D = This gives us the non-degenerate potential 



[E16] V = -=L= fa ■ ^ ■ ^ 



V^M-y^ \ X + ^/x"^ + y"^ X — yjx"^ + y"^ ) 
Here the extra constant of the motion has the form 

A2 = Mpy + d2{x,y). 

If AB = 0, by letting y —y if necessary, we can normalize so that A = 
1,-8 = and translate to get C = D = 0. This produces the non-degenerate 
potential 

[E17] r = ^ + 4 + ^ 



ZZ 2 Z\/ zz 

Here the extra constant of the motion has the form 

A2 = Mp+ + d2{x,y). 

li A = B = 0, the various possibihties have already been included under 
previous cases. 

There is one special case where an extra constant of the motion exists. If 
Ai is M^, so that di{x,y) = and M is a constant of the motion, then the 
only additional potential is 
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The constants of the motion can be taken as 

ol V ax 
A2 = Mp,-- J—^ . A, = Mpy + - ^ , X = M. 

I ^x^ + 2 ^x^ + 

The Poisson bracket relations are 

{X, A^} = -As , {X, As} = A, , {A,, A3} = XA^ , 
and these constants satisfy 

Al + Al-X'Ao-^ = 0. 
This is the well known Coulomb problem in two dimensions. 
For orbits of the fourth type, the constant of the motion Ai has the form 

A^ = Mp++p^^ + di{x,y), (14) 

corresponding to semi-hyperbolic coordinates (see Appendix ^). However, 
the only superintegrable potentials associated with this constant of the mo- 
tion have already been considered, [E4,E10,E11,E14]. 

For orbits of the fifth type, the first constant of the motion has the form 

Ai = M^+pl + d,{x,y), (15) 

corresponding to hyperbolic coordinates (see Appendix ^ . The second con- 
stant of the motion can be written in the form 

A2 = AMp+ + BMp_ + Cp\ + Dpi_ + d2{x, y). (16) 

If AB 7^ there are no cases with non-constant potential. If AB = 0, \A\ + 
> there are two new non-degenerate cases. 



v/i^ ^z{z + 2) ^z{z - 2) 
where the additional constant is 

A2 = Mp_ + d2{x,y). 
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The remaining possibihties, have been hsted earher. 

For orbits of the sixth type, the first constant of the motion has the form 

A = Mp, + di(x,y), (17) 

corresponding to parabohc coordinates (see Appendix 0). There is only one 
(non-degenerate) case that is not already listed above 0: 

[E20] V = (^a + [3^x + ,fi^+f + 7^^- yx^T^^ 

where the extra constant of the motion is 



A2 = Mpy + d2{x,y). 

Note that although this potential only separates in parabolic coordinate 
systems, it separates in more than one such coordinate system and 
hence is multiseparable. 

For orbits of the seventh type, the first constant of the motion has the form 

A^ = M^ + c'pl + di{x,y), (18) 

corresponding to elliptic coordinates (see Appendix ^), however, all superin- 
tegrable potentials separating in an elliptic coordinate system have already 
been listed. 

The last orbit on our list of equivalence classes has a typical representative 
Ai = Mp^ + di{x,y). The second constant of the motion A2 must lie on 
the equivalence class of one of the eight canonical types (^. Therefore, by 
an Euclidean group motion (including reflections), and by adding multiples 
of Aq if necessary, we can assume that the leading terms of A2 are equal to 
one of the eight representatives (^. Under this transformation Ai will be 
mapped to a constant of motion of the form Ai = Mp± + ap\ + di{x, y). For 
seven of these representatives we have already listed all possible potentials 
above. Therefore the only new case we need consider is when A2 transforms 
to A2 = Mp+ + d2{x, y). Since Ai, A2 are functionally independent constants 
of the motion, we must have either Ai = Mp_ + ap"^ + di {x, y) or a ^ 
0. Consequently the potential under consideration must admit a quadratic 
constant of the form + d^lx, y) or one that can be further transformed to 
Mpx+d^i^x, y). However, we have already listed all superintegrable potentials 
admitting a constant of one these forms. Thus there are no new potentials 
corresponding to this orbit. 

This completes the list of possible potentials involved in our problem. As 
a consequence we see that the list of 20 potentials that we have calculated 
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completely solves the problem in two dimensions of when a potential added to 
a flat space admits more than one quadratic constant of the motion. All other 
cases are equivalent to these via proper complex Euclidean transformations 
and reflections. 

3 Superintegrability on the complex two-sphere 

We can also solve the similar problem on the complex sphere. Our ba- 
sic problem is to find the superintegrable potentials for the solution of the 
Hamilton- Jacobi equation on the complex two sphere, 

H = J',+ J2 + Jl + V{x,y,z) = E, (19) 

with x,y,z subject to the constraint + y'^ + = 1, and Ji = yp^ — 
zpy^ J2 = — xpz, J3 = xpy — ypx- There are five inequivalent separable 
coordinate systems for the zero potential equation (|19|) and five different 
quadratic orbits, see Appendix 0. Typical representatives of these orbit 
classes are 

(Jl-^J2)^ 73(^1-^^2), {Ji + iJ2)-Jl (20) 
Jl Jl + r\Jl{r^±lM<l). 

We can proceed as we have in the case of the complex Euclidean plane. 
We consider one of our two quadratic constants to correspond to one of the 
representatives, (^Of ) , hence coming from a separable coordinate system in 
standard form. The potential must then have an explicit separable form in 
the appropriate coordinates. We then ask when does there exist an additional 
quadratic constant and what conditions does this impose on our potential. 
Potentials are considered as equivalent if they are related by an action of 
the complex orthogonal group 0(3), including reflections. For background 
information about this problem, see [|l^, |^. 

Unlike Euclidean space, all superintegrable potentials are multisepara- 
ble. A method for determining the type of separating coordinates from a 
given constant is described in Appendix |B| and various possibilities for each 
potential found below are listed. 

We consider first those systems that separate in horospherical coordinates. 
Thus, there is a quadratic constant of the form 

Ai = Jj + di{x,y,z). 

(We adopt the notation J± = Ji ± iJ2 and w = x + iy, w = x — iy.) In 
terms of horospherical coordinates u, v this means that the potential can be 
represented in the form 

V = f{v) + v^h{u) 
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for some functions /, h. We now assume that there is a second constant of 
the motion. It can be taken in the form 

A2 = AJ+^ + 5J3 J_ + CJ3J+ + DJ^ + v). (21) 

One can show that the case A ^ does not admit any nonconstant potentials. 
Similarly, the case A = 0,C doesn't occur, li A = C = and D 0, 
then via a symmetry transformation exp(aJ_), for suitable a we leave Ai 
unchanged and map A2 to 

A2 = DJl + Ja- 

Thus there are only two cases: 1)A = C = D = 0,3 = 1 and 2) A = B = 
C = 0,D = 1. 



a [3z 7(1 — 4z 



2^ 



[SI] v=— + '—+ ^ . 

The extra constant has the form [0 

A2 = J3J- + d2{x,y,z) . 

[S21 v = 4 + 4 + 2|. 

The extra constant has the form [0 

A2 = Jl + d2{x,y,z) . 

There is a special case of [S2] that admits an extra symmetry. 

[S3] V=^^. 

The two extra constants are of the form 

A2 = {Jl + iJ2f + d2{x,y,z), A3 = J3. 
For convenience, we will adopt a modified basis given by 



^1 = Jl + — — 7r-2 ' ^2 = J1J2 - «^ , X = J3 . 

The Poisson bracket relations are 

{X, ii} = -2^2 , {X, A2} = -Ao + X' + 2A, , 

{AuA2} = -X(2ii + a), 
with functional relation 

2 

ii(Ao - ii - X^) - - + Ao) + ^ = 0. 
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We now consider degenerate elliptical coordinates of type 2. The defining 
constant of the motion has the form 

Ai = J^J. + di{x,y,z). 

There is only one (non-degenerate) new system [0 



with constant of the motion 

^2 = ^1 + d2{x,y,z) . 

There are two special cases of [S4] that admit an extra constant of the 
motion: 

[S5] V=-^. 

where the extra constants can be taken as 

Ai = JiJ_ -, A2 = Jl + a—, X = J_. 

w w 

The Poisson bracket relations take the form 

{X,A^} = tX^-ta, {X,A2} = 2^A^, {A^, A^} = 2%X A^ 
with the functional relation 

^2-^2X^ + 0(^2-^0) = 0. 

[S6] V 



A suitable choice of basis is given by 

a X A -J J " y x-T 

The Poisson bracket relations are 

{X, A2} = -A3 , {X, A3} = A2 , {A2, As} = X{Ao - 2X2) ^ 
with the functional relation 

Al + Al + X^- AqX^ = . 
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For degenerate elliptical coordinates of type 1 the constant describing this 
system has the form 

A, = J+'-Jl + d,{x,y,z). (22) 
Two new (non-degenerate) potentials arise: 

[S7] V = , ,^ , , + „ , , , + ^. 

+ -2 z'^^Jy^ + z-^ 
The extra constant has the form 

^2 = Jl + d2ix,y,z), 

see 0. 

[S8] V" = + /^(^^ - ^) ^ 7(w + ^) 

VlP~+^ ^w{z — zy) \Jw{z + iy) 

with the second constant given by 

^2 = -^3-^1 + d2{x,y,z). 

There are no special potentials that give a third constant of the motion in 
this case. 

We now consider spherical coordinates. Here, the first constant has the 
form 

Ai = Jl + di{x,y,z). 
There is one (non-degenerate) new system |]^. The potential is 

[S9] ^=4+4+^- 

X"^ 

The extra constant has the form 

A2 = JI + d2{x,y,z). 

All of the elliptical coordinate cases have already been covered in the cases 
above. This completes the list of possible superintegrable potentials on the 
complex two-sphere. 
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4 Conclusions 



In this paper we have, in complete generahty, enumerated all potentials on 
two-dimensional complex constant curvature spaces for which there is more 
than one constant of the motion that is quadratic in the momenta. For 
each pair of constants of the motion, whose leading terms are second order 
in the enveloping algebra of the Lie symmetry algebra of the free particle 
Hamiltonian, we find a pair of coupled second order linear partial differential 
equations satisfied by the potential function. The key to making our approach 
practical is that when one of the constants of the motion corresponds to a 
separable coordinate system, we can explicitly (and simply) solve one of these 
PDEs in this coordinate system, and merely have to substitute the solution 
into the second equation. 

One can see by inspection of Tables ^ and ^ that each of these cases 
(except one) is multiseparable, i.e., separation is possible in at least two 
coordinate systems. The one counterexample in flat space ([E15]) still sepa- 
rates in one system. These tables also show that each potential listed can be 
uniquely identified by its list of associated equivalence classes of quadratic 
constants. This serves to confirm that they are indeed distinct potentials, 
unrelated by group motions. 

We also observe that whenever there is more than one extra quadratic 
constant, a first order constant can be found. Further, the non-degenerate 
potentials found in 0] that are not related to a degenerate potential by 
group motions are those for which the additional constants are genuinely 
second order, i.e. no first order constant exists. 

Note further that for a non-degenerate potential in fiat space we can pre- 
scribe V, Vx, Vy, Vyy arbitrarily at any regular point (a;o,yo) and these values 
determine V"(x, y) uniquely. These potentials correspond to exactly three 
functionally and linearly independent constants of the motion. For a degen- 
erate potential with an extra (linearly independent) constant of the motion 
the additional constant implies a relationship between Vx,Vy at any regular 
point; hence that all first, second and higher order derivatives of V{x,y) can 
be expressed in terms of a single first derivative, say V^. Thus for all these 
potentials we can prescribe V, arbitrarily at any regular point and these 
values determine V{x, y) uniquely. It follows that except for the exceptional 
case ([E15]) the superintegrable potentials depend on exactly 4 or 2 parame- 
ters. Analogous comments hold for the complex two-sphere, except that here 
there is no exceptional case. 

What is exceptional about ([El 5])? This is the only case where one 
cannot solve for V^x — Vyy and V^y as linear combinations of Vx,Vy. Thus 
the potential must be degenerate. Indeed this potential, although it depends 
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on an infinite number of parameters, must have the form V{z) so Vx and Vy 
cannot be prescribed independently at a point. Furthermore, the potential 
is not uniquely determined by the values of V, Vx,Vy, Vyy at a point. 

We give in this paper, and preceding papers, the structure of the classical 
quadratic algebras in almost all cases. We intend to perform a comprehensive 
study of the corresponding quantum algebras associated with the Schrodinger 
equation at a later date. 



Each coordinate system in which the Hamilton- Jacobi equation is separable 
on E2fi is characterized by a constant quadratic in the momenta. Coordinate 
systems that are related by Euchdean group motions belong to the same 
family and hence a given family of coordinates (e.g. polar coordinates) is 
associated with an equivalence class of quadratic elements in the enveloping 
algebra of e(2, C). Two elements are equivalent if one can be transformed 
into to other by a combination of scalar multiplication, addition of multiples 
of + and Euclidean motions (including reflections) . One equivalence 
class (listed below) is not associated with a separating coordinate system. 
The following can be taken as a representative list of coordinate systems and 
corresponding constants. 

1. Cartesian coordinates. 



A Separable coordinates in ^2,c 




2. Light cone coordinates. 



z — X -\-iy , z — X 



3. Polar coordinates. 



xs — r cos ^ , y5 = r sin ^ , 



4. Semi-Hyperbolic coordinates. 



xsH ^i{w - uf + 2i{w + u) , ysH 



{w - uf + 2{w + u) , 



L = M{px + ipy) + {px - iPyf ■ 
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5. Hyperbolic coordinates. 

r^ + s^ + r^s^ + - r^s^ ,^2 / \2 

" 2^s ' " * 2^^^ ' ^ ^ ^ + (Px + ■ 

6. Parabolic coordinates. 

xp = iv, yp = ^(C^ - , L = Mp^ . 

7. Elliptic coordinates. 

xe = c^{u — l){v — 1) , = cV—uv , L = + c^p^ . 

8. No separation. 

No corresponding separable coordinates, L — M{px + ipy) . 

The following facts are useful in determining to which class given constant 
belongs. 

• Translations leave p^ and Py unchanged and for any A and B a trans- 
lation can be found that has the effect 

M + Ap^ + Bpy. 

• Rotations leave M fixed and one can be found that for any A and B 
has the effect 

Ap^ + Bpy one of p^, p+ or p_. 

• A rotation can be found that for any Ay^O has the effect 

p.^Xp,. and p^^ip^. 

• The reflection y ^ —y has the effect 

M -M , p+ < — > p_ . 

For each superintegrable potential [El— 20], all linear combinations of the 
given quadratic constants must be considered in order to determine which 
equivalence classes are represented, and hence in which families of coordinates 
systems it will separate. 
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For example, the potential [E18] has constants with leading part L = 
AM"^ + BMpx + CMpy + Dlp^+Py). From this we can see immediately that 
polar and parabolic coordinates will separate this Hamiltonian, and that the 
non-separating constant Mp+ can be generated. Further, M^+2Mpy + 2(p^ + 
Py) +pI under a translation that maps M — > M — py and hence [E18] 

separates in an elliptic coordinate system. Lastly, there exists a translation 
mapping + 2iMp^ — ^ + p"^ , and hence the Hamiltonian separates in 
hyperbolic coordinates. 

The results of similar reasoning for the other potentials are summarized 
in table 0. 



B Separable coordinates in 82,0 

As for Euclidean space, coordinates separating the Hamilton- Jacobi equation 
on the two-sphere correspond to constants that are quadratic in the elements 
of the Lie algebra of its symmetry group 0(3). Coordinates belong to the 
same family if one can be transformed to the other by a rotation or reflec- 
tion. On the complex two-sphere, unlike complex Euclidean space, every 
quadratic constant, other than a multiple of the Hamiltonian, corresponds 
to a separating coordinate system. 

The separable coordinates on the complex two-sphere and their charac- 
terizing constants are: 

1. Spherical coordinates. 

X = sin 9 cos (p , y = sin 9 simp , 

z = cos 9 , L = J? 



3 



2. Horospherical coordinates. 



X 





V 

3. Elliptic coordinates. 



z =— , L = ( Ji - iJ2) 



2 (ru — l){rv — 1) 2 r{u — l){v — l) 
1 — r r — 1 

z"^ = ruv , L = + r J2 . 
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no. of distinct eigenvalues 


dimker((/)) 


Spherical 


2 


1 


Horospherical 


1 


1 


Elliptic 


3 





Degenerate elliptic type 1 


2 





Degenerate elliptic type 2 


1 






Table 1: Invariants used to identify coordinate systems on 82^0 



4. Degenerate Elliptic coordinates of type 1. 

{u^ + l){v^ + l)' L-[J, + iJ,) -J3. 

5. Degenerate Elliptic coordinates of type 2. 

l(M2 + t;2)2 



X + ly = —luv , X — ly 



2 uv ' 



4 vP'v^ ' 



The action of the symmetry group on a general quadratic constant is 
not as easily described as for -£'2,0 • To determine the equivalence class to 
which a given quadratic element L belongs it is more convenient to note that 
the number of distinct eigenvalues of L, as a quadratic form in the Jj, and 
the dimension of the kernel of the map on first order elements, cj) : X ^ 
{X, L}, are both invariant under group motions and addition of multiples of 
the Casimir Jl + + J^. Table [1| gives the correspondence between these 
invariants and families of coordinate systems on 82,0 ■ 

Just as for E2^Ci by considing a general linear combination of constants for 
each potential [SI— S9], the corresponding families of separable coordinates 
can be determined. The results are summarized in table ^ 
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El 


E2 


E3 


E4 


E5 


E6 


E7 


E8 


E9 


ElO 


Cartesian 


X 


X 


X 


X 


X 


X 










Light Cone 






X 


X 


X 




X 


X 


X 


X 


Polar 


X 




X 






X 




X 






Semi- Hyperbolic 








X 












X 


Hyperbolic 






X 








X 


X 






Parabolic 




X 






X 


X 






X 




Elliptic 


X 




X 






X 


X 








Non-Separating 








X 
















Ell 


E12 


E13 


E14 


E15 


E16 


E17 


E18 


E19 


E20 


Cartesian 






















Light Cone 


X 


X 


X 


X 


X 












Polar 








X 




X 


X 


X 






Semi- Hyperbolic 


X 




X 
















Hyperbolic 




X 




X 






X 


X 


X 




Parabolic 






X 






X 




X 




X X 


Elliptic 




X 








X 




X 


X 




Non-Separating 


X 


X 


X 


X 


X 




X 


X 


X 


X 



Table 2: Separating coordinate systems for superintegrable potentials in com- 
plex two-dimensional Euclidean space. Potentials possessing a quadratic con- 
stant equivalent to Mp+ are indicated in the line labelled 'Non-Separating'. 





SI 


S2 


S3 


S4 


S5 


S6 


S7 


S8 


S9 


Spherical 




X 


X 


X 


X 


X 


X 




X 


Horospherical 


X 


X 


X 




X 










Elliptic 


X 




X 




X 


X 


X 


X 


X 


Degenerate elliptic 1 




X 


X 


X 


X 


X 


X 


X 




Degenerate elliptic 2 


X 






X 


X 


X 









Table 3: Separating coordinate systems for superintegrable potentials on the 
two-dimensional complex sphere. 
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